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Abstract: 

Parton percolation provides geometric deconfmement in the pre-equilibrium stage of 
nuclear collisions. The resulting parton condensate can lead to charmonium suppression. 
We formulate a local percolation condition viable for non-uniform collision environments 
and show that it correctly reproduces the suppression observed for S — U and Pb — Pb 
collisions at the SPS. Using this formulation, we then determine the behavior of J/ip sup- 
pression for In — In collisions at the SPS and for Au — Au collisions at RHIC. 
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1. Introduction 



In recent years, the role of the partonic initial stages of high energy nuclear collisions 
has attracted increasing interest. While the ultimate aim of the experimental program is 
the production of the quark-gluon plasma predicted by statistical QCD, it is becoming 
more and more evident that a thermalized medium consisting of quarks and gluons can be 
produced only if already the initial state of the collision provides the conditions necessary 
for deconfmement and thermalization. Dilute initial state parton configurations do not 
lead to the colour connection needed so that partons from different collisions can combine 
to form a collective medium. On the other hand, if the initial state of primary partons 
is sufficiently dense, cluster percolation will set in, leading to a condensate of connected 
and hence interacting partons, which are no longer associated to any "parent" hadrons. 
Thus parton percolation is a geometric, pre-equilibrium form of deconfmement PQI2]; it 
must occur in any description based on partons of an intrinsic transverse momentum. 

Subsequent interactions can thermalize this interconnected partonic system, turning 
it into a quark-gluon plasma, so that parton percolation constitutes an essential prere- 
quisite for QGP formation. Moreover, pre-equilibrium deconfmement leads to a further 
important question. Much of the investigation of high energy nuclear collisions is devoted 
to the search for quark-gluon plasma signatures. However, some of the features observed 
in nuclear interactions might be determined by the system present before any equilibration 
has occurred, and they could be independent of a subsequent thermalization. We must 
therefore investigate observable consequences of parton percolation, and if possible, show 
how they can be distinguished from those due to QGP formation. 

One of the features of particular interest in this connection is the suppression of 
J I ip production in nuclear collisions, predicted as deconfmement signal for thermal media 
quite some time ago [Sj. If the parton condensate formed through percolation contains 
partons hard enough to resolve the produced charmonium states, these partons can also 
dissociate charmonia, so that the onset of J/ip suppression could coincide with that of 
parton percolation. It was shown in a first study jl] that for Pb — Pb collisions at the 
CERN-SPS as function of centrality the resulting thresholds appear quite reasonable. A 
more quantitative comparison, which can then also provide the basis for further predic- 
tions, requires a formulation of percolation in a non-uniform environment, as given by the 
partonic source profile in nuclear collisions, and adapted to the local nature of charmonium 
probes. The aim of this paper is first to formulate percolation for such conditions, then to 
compare the results to existing SPS data (S — U and Pb — Pb), and finally to present the 
predictions of this approach for the forthcoming J/ijj production experiments at CERN 
and BNL. All results are obtained by Monte Carlo simulations taking into account the 
centrality dependence of the collisions. In closing, we shall consider some features which 
could distinguish J/ ip suppression by parton percolation from a suppression in a thermal 
medium. 

2. Local Parton Percolation Conditions 

Consider a flat two-dimensional circular surface of radius R (the transverse nuclear 
area) , on which N small discs of radius r<i? (the transverse partonic size) are randomly 
distributed, allowing overlap. With increasing density n = N/nR 2 , clusters of increasing 
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size appear. The crucial feature is that this cluster formation shows critical behaviour: in 
the limit N — > oo and R — > oo, with n finite, the cluster size diverges at a certain critical 
density n = n c . This percolation threshold is given by 

n c= — o C 1 ) 

the critical value z/ c ~ 1.13 of the 'filling factor' v = n(r/R) 2 is determined by numerical 
studies. For finite N and R, percolation sets in when the largest cluster spans the entire 
circular surface from the center to the edge. The resulting cluster growth is illustrated 
in Fig. for a ratio r/R = 1/100; we show the percolation probability, defined as the 
relative size of the largest cluster, as function of the filling factor v. Because of overlap, a 
considerable fraction of the surface is still empty at the percolation point; in fact, at that 
threshold, only 1 — exp{— u c } ~ 2/3 of the surface is covered by discs. 




Figure 1: Cluster growth and corresponding derivative as function of the overall filling 
factor v (a) and of the cluster filling factor r\ (b). 



In high energy nuclear collisions, the incident partons are distributed on a plane trans- 
verse to the collision axis. However, since they originate from the nucleons within the 
colliding nuclei, this distribution is highly non-uniform, with more nucleons and hence 
more partons in the center than towards the edge of the transverse nuclear plane. If the 
basic surface is not flat, it is still possible to define the percolation threshold as the point 
at which the cluster size shows singular behavior in the limit of large R and N [2] . Such a 
global definition is, however, not always the most useful. Hard probes, such as quarkonia, 
probe the medium locally and thus test only if it has reached the percolation point and 
the resulting geometric deconfinement at their location; they cannot register the global 
features of the entire cluster area. It is thus necessary to define a more local percolation 
criterium, and this is in fact quite straight-forward. 

As mentioned, at the percolation point, exp{— u c } ~ 1/3 of the surface remains empty. 
Hence the disc density in the percolating cluster must be greater than (1.5 v c /-Kr 2 ). Nu- 
merical studies show that in fact percolation sets in when the density m of constituents 
in the largest cluster reaches the critical value 

m c = \ (2) 
iir 2 



3 



with T] c ~ 1.72, slightly larger than 1.5 v c . This result provides the required local test: if 
the parton density at a certain point in the transverse nuclear collision plane has reached 
this level, the medium there belongs to a percolating cluster and hence to a deconfmed 
parton condensate. In Fig. QJd, we illustrate how the percolating cluster grows as function 
of the cluster filling factor 77 = m(r/R) 2 , again with r/R = 1/100. 

Let us recapitulate: there are two equivalent criteria which can be used to specify 
the onset of percolation. In the 'global' definition, it occurs when the total number 
of partons distributed over the entire transverse nuclear area reaches 1.13/irr 2 . In the 
'local' definition, percolation sets in when the parton density in the largest cluster reaches 
1.72/nr 2 . Which of the two is preferable in a given case depends on the physics question 
addressed, and if we want to study local phenomena such as J / ip suppression, it is clearly 
the local parton density which is relevant, not the average over the entire transverse area. 

We now turn to the implications of the approach to nuclear collisions j3j; for illus- 
tration, we concentrate for the moment on central A — A interactions at a c.m.s. energy 
y/s per nucleon-nucleon collision. The distribution of nucleons in the colliding nuclei is 
specified by a Glauber calculation using Woods-Saxon nuclear distributions 5J; this pro- 
vides the density n s (A) of nucleons in the transverse collision plane. The parton content 
of a nucleon is given by parton distribution functions dN q (x,Q 2 )/dy determined in deep 
inelastic scattering experiments; here x denotes the fraction of the nucleon momentum 
carried by the parton and Q the momentum resolution scale. At central c.m.s. rapidi- 
ty y — 0, we have x = kx/y/s, where kx denotes the average transverse momentum of 
the parton. In nuclear collisions, kx defines the transverse size of the partons and thus 
also sets the resolution scale, kx — Q- Using these quantities, we have for the 'global' 
percolation condition in nuclear collisions 



it determines for what value of A at a given y/s percolation sets in, and it specifies the value 
Q = Q C (A, y/s) at the onset point. Using the nuclear source density n a (A) together with 
the parton distribution function, we can also calculate the density of partons m(A, Q, y/s) 
for the largest cluster in the transverse collision plane. The relation 



then provides the 'local' parton percolation condition to specify the onset values of A, y/s 
and Q. These considerations can be extended to non-central collisions in a straight- 
forward fashion |3J |S] . 

It is not a priori evident that the threshold values v c and r/ c for the non-uniform 
distribution specified by the Woods-Saxon distributions remain the same as for a uniform 
distribution, nor is it clear how large finite size corrections are at given values of r/R. 
We have therefore calculated the pseudo-critical values (derivative peak positions) of v 
and 77 for R/r = 20, 50 and 100, using the Woods-Saxon distribution for central Pb — Pb 
collisions to determine the parton density as well as the average overall interaction region. 
Finite-size scaling techniques from statistical physics then provide the asymptotic limit; 





m c (A, Q c , y/s) 




(4) 
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in this way, we have verified that v c and r\ c for this non-uniform distribution in fact agree 
with the above mentioned values for the uniform case. Moreover, we found that for the 
mentioned values of R/r, the deviations from the asymptotic critical values are less than 
2%. 

2. Comparison to Pb-Pb and S-U Data 

We now want to study parton percolation for CERN-SPS conditions and compare the 
resulting consequences for J /if} suppression to the data of the NA38/NA50 collaborations. 

The parton content of a nucleon is obtained from GRV94 parton distribution functions 
jHj, since these are available for the required kinematic range; for more details, see [I]. 
With this specified, we calculate the size of the largest cluster as well as the overall size 
of the interaction region for different centralities (the latter is here defined by the ratio 
of the number to the density of wounded nucleons). For Pb — Pb collisions at ^/s = 11 A 
GeV, we show in Fig. |2K the relative size of the largest cluster as function of the filling 
factor 77. It is seen that the cluster size grows rapidly in a narrow band around rj c , and at 
the percolation onset, the percolating cluster covers about 1/3 of the total collision area. 
In Fig. |2t>, we show the cluster filling factor 77 as function of centrality, measured by the 
number of participating nucleons. The critical value for parton percolation, rj c = 1.72, is 
found to be attained for N part ~ 125. 
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Figure 2: Percolation onset in Pb-Pb collisions at the SPS. 



The parton distribution functions specify both the number of partons and the resolu- 
tion scale Q c at the percolation point. For Pb—Pb collisions at \fs = 17.4 GeV, Q c ~ 0.7 
GeV. The scales of the charmonium states Xc and if)' , as determined by the inverse of 
their radii calculated in potential theory, are around 0.6 GeV and 0.5 GeV, respectively. 
Since the parton condensate can thus resolve these states, we assume, following 0j, that it 
will dissociate them. To be specific: we assume that all \ c and if> f states formed inside the 
percolating cluster disappear; the location of their formation is determined by the collision 
density as given by the Glauber formulation [3J |3] . The first onset of J / if> suppression 
in Pb — Pb collisions at the SPS should therefore occur at iVp art ~ 125, where the J/^'s 
due to feed-down from Xc and ip' states in the percolating cluster are eliminated. Directly 
produced J/ip's survive because of their smaller radii (leading to a scale of 0.9 - 1.0 GeV), 
as do those coming from excited states outside the percolating cluster. We emphasize that 
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with parton distribution functions and charmonium radii given, the description contains 
no adjustable parameters. 

As noted, the parton condensate formed at the onset of percolation has a resolution 
scale Q c ~ 0.7 GeV and thus cannot resolve the ground state charmonium J/ifj with its 
smaller radius. The dissociation of directly produced J/ip's thus requires more central 
collisions, which lead to a higher parton density and simultaneously to a better resolution, 
i.e., to an increase of Q c . The precise onset point is here more difficult to determine, since 
it is fixed by the value of the J/ip radius. For rj/^ = 0.20 fm, we have Q ~ 1.0 GeV 
and obtain the centrality dependence shown in Fig. El indicating that the suppression of 
directly produced J/ijj's starts at iVp art ~ 320. However, here a note of warning has to 
be added. The Xc and ip' states can be resolved by the parton condensate when it is first 
formed; the values of their intrinsic scales do not enter explicitly, since they are below 
the resolution scale Q c of the condensate. The onset point for the suppression of the IS" 
ground state is directly determined by the radius of that state, and since the dependence 
of parton density on centrality (see Fig. EJ) is quite flat, small changes in rj/^ lead to large 
changes in the threshold value. For instance, rjj^ = 0.22 fm (Q ~ 0.9 GeV) shifts the 
onset from A^ part ~ 320 to about 190, as also shown in Fig. El 
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Figure 3: Percolation onset for directly produced J/^'s in Pb-Pb collisions at the SPS. 

We now combine both thresholds and use the nuclear collision profile to determine the 
distribution of charmonium production in the transverse plane 0, in order to determine 
which are formed in the cluster and which are not. As function of centrality, this leads to 
the survival pattern shown in Fig. |3J 

It is clear that various theoretical and experimental sources (fluctuations, resolution, 
binning, etc.) will lead to some smearing of the theoretical survival pattern. To obtain 
some idea of what this may lead to, we have folded our calculated distribution .S^part) 
with a Gaussian smearing function, 

S(N part ) = I dN part S(iV part )-=L- exp{-(iV part - A part ) 2 /2a 2 }, (5) 

using a relative resolution of cr/N paxt = 15 % at the respective thresholds. The result is 
also shown in Fig. 0J 

The centrality dependence of the J / ip survival has been measured in Pb-Pb collisions 
at the SPS using two independent variables, the associated transverse energy, Er, and 
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Figure 4: J /if) survival pattern in Pb-Pb collisions at the SPS, as function of centrality. 



the energy of the undeflected nucleons, -Ezdc- The first effectively measures the number 
of hadronic secondaries in the transverse plane, the second the number of participating 
nucleons. We want to study anomalous charmonium suppression, i.e., effects above and 
beyond the normal suppression already observed in p-A collisions, due to absorption of 
the nascent charmonium state in normal nuclear matter jHJ. To determine this normal 
suppression, we use a recent analysis [H] ofp—A and S — U data by the NA50 group, 
which finds a somewhat weaker nuclear absorption, ad — 4.3 ± 0.6, than previous studies. 
Normalizing the Pb-Pb data to the normal nuclear suppression thus specified leads to the 
points shown in Fig. where they are compared to our parton percolation suppression 
pattern. We see good agreement between the calculated curve and the data points. We 
recall, however, that the second threshold value depends crucially on the J /if) radius. 
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Figure 5: J /if) survival pattern in Pb-Pb collisions at the SPS, as function of centrality, 
determined from Et (left) and from Ezdc (right). 



Next we want to check the behaviour predicted for S—U collisions at the SPS, since this 
shows at best only very slight anomalous suppression 1 . In Fig- El we show the percolation 

initially, it was claimed that the data are fully accounted for by normal nuclear absorption. However, 
the lower absorption cross sections obtained in the mentioned new p-A studies allow a slight anomalous 
suppression for the most central S—U collisions. 
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behaviour, and in Fig. [7[ the resulting survival pattern together with the data [TT]. The 
threshold obtained here corresponds to the general onset of percolation; the resolution 
scale of the direct J ftp cannot be reached in S— U collisions. The corresponding smeared 
form of the calculated distribution is also included. 
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Figure 6: Percolation onset in S-U collisions at the SPS. 
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Figure 7: J ftp survival pattern in S-U collisions at the SPS, as function of centrality. 



We close this section with some general comments on the suppression of charmonium 
production. We have here considered parton percolation as mechanism for anomalous 
J ftp suppression, defined as the reduction of the observed production rate beyond that 
predicted by normal absorption in standard nuclear matter. In doing so, we neglect 
other conceivable suppression mechanisms, in particular dissociation by (non-percolating) 
secondary hadronic or partonic comovers. This is not expected to play a significant role 
for the suppression of IS states, i.e., directly produced J/ip's: the large binding energy 
(650 MeV) makes a dissociation by confined partons very unlikely [T2]. However, since 
about 40 % of the observed J ftp's come from the decay of higher excited states (feed-down 
from Xc and ip' production), we have to consider the fate of these as well. The Xc (~ 250 
MeV) is still rather strongly bound, but the much smaller ip' binding energy (~ 60 MeV) 
means that any interaction will cause break-up. It is thus not surprising that the ip' suffer 
an anomalous suppression even in the most peripheral S-U interactions • This is most 
likely not related to any parton deconfinement and has to be addressed separately. 
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3. Predictions for NA60 and RHIC 



We now turn to the forthcoming In— In collisions at the SPS. The percolation pattern 
is shown in Fig.|S] the general percolation threshold is reached at N pait ~ 140. In contrast 
to Pb — Pb collisions, the threshold for directly produced J/ip's here remains above the 
produced density even for the most central collisions (but recall the caveat concerning its 
dependence on the J /if) radius). The resulting survival pattern for In — In collisions at 
the SPS is given in Fig. |Hl including also the smeared form. We see that the most central 
30 % of the collisions show a suppression of the J/if)'s produced by feed-down from higher 
excited charmonium states. 
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Figure 8: Percolation onset in In-In collisions at the SPS. 
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Figure 9: J /if) survival pattern in In-In collisions at the SPS, as function of centrality. 



Finally we address the suppression pattern for Au — Au collisions at RHIC, with the 
percolation pattern shown in Fig. Here the increased parton density shifts the onset 
of percolation to a higher resolution scale, so that from the threshold on, all charmonium 
states are suppressed. This leads to a single step suppression pattern, starting at N part ~ 
90; it is given in Fig. ^2 Recall that we are here specifying the predicted anomalous 
suppression, beyond the normal suppression resulting from (pre-resonance) absorption in 
standard nuclear matter. For various reasons, this normal suppression may depend on 
the collision energy and can thus be quite different from that observed at SPS energy. 
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To determine the anomalous suppression at RHIC, it is therefore crucial to first measure 
the normal absorption in p— A (or d— A) collisions. The threshold for Au — Au collisions 
as obtained here occurs for considerably more central collisions than found in previous 
studies j2], in which the dependence of parton distribution functions on the resolution 
scale had not taken into account. 
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Figure 10: Percolation onset in Au-Au collisions at RHIC. 
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Figure 11: J/ip survival pattern in Au-Au collisions at RHIC, as function of centrality. 



In summary: if the origin of the two-step suppression pattern observed in the Pb—Pb 
SPS collisions is indeed the onset of parton percolation at different resolution scales, as 
discussed, then we expect a single step pattern both in In — In collisions at the SPS and 
in Au — Au collisions at RHIC. In the former, the step is due to the suppression of the 
higher excited states only, with directly produced J/^'s surviving; in the latter, there is 
a common threshold for all charmonium states and hence only one step. 

4. Discussion 

We have presented the J/ip suppression pattern as obtained from parton percolation 
under the assumption of complete charmonium dissociation within a percolating cluster of 
sufficient resolution scale, full survival in media below percolation density or of insufficient 
resolution scale. In this closing section, we want to summarize the uncertainties of such 
an approach and compare its outcome to that of a thermal description. 
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The onset of parton percolation is specified directly through the parton distribution 
functions (PDF's) as obtained from an analysis of deep inelastic lepton-hadron data. Un- 
certainties in the PDF's produce some uncertainty in the precise values of the percolation 
threshold, but not in the onset as such. As mentioned, percolation must occur in any ap- 
proach based on partons with an intrinsic transverse momentum. The assumption tested 
by comparing our results to J / ip data is that for charmonium suppression in nuclear colli- 
sions, the primary parton density is the crucial variable and the formation of a percolating 
cluster the relevant mechanism. 

As noted, the onset of the suppression of directly produced J/^'s at SPS energy is quite 
sensitive to the choice of charmonium parameters (and on that of the parton distribution 
functions, for which slight changes would also produce considerable threshold shifts). On 
the other hand, the onset of the first step in SPS Pb—Pb pattern, the predicted single step 
in SPS In— In collisions, and the predicted single step in RHIC Au—Au collisions should be 
quite robust. This should help in identifying the origin of the second step in SPS Pb—Pb 
collisions. It has been suggested [Hj that this drop is largely due to transverse energy 
fluctuations for the most central collisions; this would imply that it should be absent in 
zero-degree calorimeter studies, where (see Fig. EJ) it seems to still occur to some extent 
[Hj. A further test should come from the SPS In— In and the RHIC Au—Au data, where 
the existence of a second step could not be attributed to parton percolation. 

Finally we turn briefly to the question of the distinguishing the initial state char- 
monium suppression presented here from a suppression occurring at a later thermalized 
stage. We had here assumed that the primary parton configuration determines the pat- 
tern of J/ip suppression. What would one expect if the initial parton configuration leads 
to subsequent thermalization, with the formation of a more or less equilibrated quark- 
gluon plasma, and it is this QGP which causes the observed J/i[) suppression? The most 
recent (quenched) lattice QCD studies [13 HE] indicate that in a hot thermal medium 
of quarks and gluons, higher excited charmonium states are dissociated at or below T c . 
In contrast, the ground state J/ip survives to considerably higher temperatures; present 
calculations still show a clear J/ip signal at least up to T = 1.5 T c . This would imply that 
thermal dissociation removes the feed-down contributions from if)' and Xc states for energy 
densities around or below 1 GeV/fm 3 , while energy densities of more than 6-8 GeV/fm 3 
are required for the melting of directly produced J/i/j's. It seems difficult to accommo- 
date the existing SPS data with these thresholds, while the pattern resulting from parton 
percolation can do so quite naturally. For Au — Au collisions at RHIC, thermal melting 
would suppress the higher state feed-down contributions for practically all centralities, 
while the directly produced J/ip's disappear at best for the most central collisions. We 
thus conclude that the charmonium suppression thresholds from initial state and from 
thermal suppression are quite different also at much higher energy. 
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